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Abstract 



In order to highlight the onset of chaos in the PuUen-Edmonds model 
a quantal analog of the resonance overlap criterion has been examined. A 
quite good agreement between analytical and numerical results is obtained. 
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Recently great interest has been shown in the so-called quantum chaos, 
or more precisely quantum chaology, i.e. the study of quantal systems which 
are chaotic in the classical limit — > [1] . 

The aim of this work is to apply a quantal analog of the resonance cri- 
terion [2] to the Pullen Edmonds model [3], in order to highlight the onset 
of chaos in a schematic model. A requisite for the occurrence of quantum 
chaos is that the spectrum behaves in an irregular way as a function of the 
coupling constant: in particular if multiple avoided crossings occur (see for 
example [4] and references quoted therein). In fact, as stressed by Berry 
[5], in quasi-integrable systems the presence of many quasi-crossings leads to 
a chaotic behaviour. Using the semiclassical quantization we calculate the 
critical value Xc of the coupling constant corresponding to the intersection 
of two neighboring quantal separatrices. This criterion is, in some sense, the 
quantal counterpart of the method of overlapping resonances developed by 
Chirikov [6]. We point out however that the Pullen Edmonds hamiltonian 
is the perturbation of an intrinsically resonant hamiltonian (see e.g. [7] and 
references quoted therein) which yields intrinsically degenerate quantum un- 
perturbed levels: hence we first have have to extend to this situation the 
argument of [2] , valid in the case of accidental degenaracy. 

As discussed by Pullen and Edmonds [4], the hamiltonian: 

H-l{pl+pl) + l{ql + ql) + Xqlql (1) 
exhibits a transition from order to chaos as a function of the coupling constant 
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X- We introduce the new variables (/, 9) by the canonical transformation: 

Qi = yflli COS Oi 

^ Pi = ^/2TiSm9i. 

Then (1) becomes: 



i = 1,2. 



(2) 



H = h + l2 + Axhh cos^ Oi cos^ 62. 



(3) 



By the new canonical transformation: 

A = h + h 
A2 = h-h 
H can be written: 



01 + 02 
01 - 02, 



(4) 



= + x(>l? - Al) COS^ (01 + 02) COS^ (01 - 02). 



(5) 



We now eliminate the dependence on the angles to order by resonant 
(or secular) canonical perturbation theory [8]. First we average on the fast 
variable (0i). This yields: 

1 r'^'" 1 

— / d0i cos^ (01 + 4)2) cos^ (0, - 02) = - (2 + cos 402) , (6) 



and: 



Ha = A^ + ^{Al- Al){2 + cos 402). (7) 



The dependence on 02 is now eliminated by a second canonical transforma- 
tion. The Hamilton-Jacobi equation for the perturbation part is indeed: 

dS 



[^?-(^)'](2 + cos402)=X, 



dS 



± 



Af(2 + cos402) -K 
2 + cos 402 



(8) 
(9) 



and thus the Hamiltonian (7) becomes: 



H = B, + ^K{B^,B2), (10) 

o 



where: 



B, = A„ B, = -fi^,-. (11) 

It appears from the structure of equation (8) that the motion of our system 
is similar that of a simple pendulum: 

< K < Bf rotational motion 
K — Bf separatrix (12) 

B\ < K < 3Bl librational motion. 

On the separatrix: 

Bl{2 + cosA(f)2)^K, (13) 

and: 



^- = ^?/''-V 2 + COS4. ■ (") 



where: 



a — h — J rotational motion 

a — (f)-{K, Bi), b — (j)+{K, Bi) librational motion 



(15) 



with: 



(f)±{K, Bi) = ±- arccos(-2 - 2). (16) 



K 

"4 'Bf 

The semiclassical quantization on the Hamiltonian (10) may be used. Set: 



Bi = mih, B2 = m2h; (17) 
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then [9], up to terms of order h, the quantum spectrum is: 

X 

Emi,m2 = rnih + -K{mih, m2h), 

o 

where K is imphcitly defined by the relation: 

2 fb 



(18) 



, 2 l(mih)^(2 + cosAx) - K 



and: 



a = -f, h = l 



0< K < {mihy 
a^(l)4K,Bi), b^(l)+{K,Bi) {mihy < K < 3{mihy 

On the separatrix, where K = {mih)'^, m2 — iami, with: 



(20) 



a — — [ ^ dx\ — 

TT i-f V 2 



+ cos Ax 



+ cos 4x 



(21) 



Figure 1 shows K vs m2. 

By (18) it is not difficult to calculate the critical value Xc of the cou- 
pling constant corresponding to the intersection of the separatrices of two 
neighboring multiplets: 

X X 
(mi + l)h + —K[{mi + l)h, a{'mi + l)h] = rriih + —K{mih, arriih), (22) 

8 8 



and so: 



Xc = 



-8h 



. (23) 

K[(mi + l)h, a{mi + l)h] — K{m-Ji, amih) 

The denominator of (23) can be evaluated by the Taylor expansion and fi- 
nally: 

(24) 



Xc 



-8 

dK dlT 
dBi " dB2 



B\ =mih,B2=am2h 
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K is implicitly defined by the relation: 



Bo- 



TT 



dx\ 



lBf{2 + cosAx) - K 
2 + cos 4x 



or: 



F{Bi,B2,K) = B2-^Bi,K) = 0. 
As a function of Xc can be written: 



0, (25) 



(26) 



8|f 

dK 



a — 



d<S> 
dBi 



(27) 



Bi =mih 



where: 



9$ 
dK 



1 rf 



da; 



TT J-f 



9$ 



^(2 + cos4x)[Sf (2 + cos 4a;) - K] 

2 ri 



dBi TT 

It is not difficult to show that 



dx, 



Bf{2 + cos4:x) 
\\ Bf{2 + cos4x) - K' 



Xc 



rriih 



(28) 



(29) 



Figure 2 shows Xc vs mih. If we take the energy E' ~ mi/i we obtain Xc — 
4/E, a quantal rule for the onset of classical chaos. 

The transition order-chaos in systems with two degrees of freedom may 
be also studied by the curvature criterion of potential energy [10,11]. It 
is however important to point out that in general the curvature criterion 
guarantees only a local instability and should therefore be combined with the 
Poincare sections [14]. For a fuller discussion of this point see [12]. 
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At low energy, the motion near the minimum of the potential V{qi, ^2) = 
+ Q2) + XQiQh where the curvature is positive, is periodic or quasiperi- 
odic and is separated from the region of instability by a line of zero curva- 
ture; if the energy is increased, the system will be for some initial conditions 
in a region of negative curvature, where the motion is chaotic. In accor- 
dance with this scenario, the energy of order— >chaos transition Ec is equal 
to the minimum value of the line of zero gaussian curvature K{qi, ^2) of 
the potential-energy surface of the system. For our potential the gaussian 
curvature vanishes at the points that satisfy the equation: 

MM - <aSs'^ = ^ '^^'^'^ ^ - = 

It is easy to show that the minimal energy on the zero-curvature line is given 
by [13]: 

E,^Vmin{K^Q,n) = ^, (31) 



and occurs at ^1 = This result is in rough agreement with equation 

(29) and with the numerical calculations of Poincare sections [14] obtained 
by authors of reference [3]. The discrepancy between (29) and (31) is due 
to the fact that the curvature criterion describes the onset of local chaos 
associated to quasi-crossings between levels of the same multiplet (see also 
reference [3]) meanwhile the quantum resonance criterion describes the on- 
set of widespread chaos associated to quasi-crossings between separatices of 
different multiplets. 



8 



The authors are greatly indebted to Dr. Stefano Isola for many useful 
discussions. 
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Figure Captions 



Figure 1: The semiclassical energy K vs 1712 for a fixed mi: < < 1 
rotational motion (a) and 1 < i^T < 3 librational motion (b). 
Figure 2: The critical value Xc vs rriih. 
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